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There is evidence of an overdensity of strong intervening Mgll absorption line systems distributed along the lines of sight towards 
GRB afterglows relative to quasar sight-lines. If this excess is real, one should also expect an overdensity of field galaxies around 
GRB sight-lines, as strong Mgll tends to trace these sources. In this work, we test this expectation by calculating the two point 
angular correlation function of galaxies within 120" (~ 470 tCj\ Kpc at (z) ~ 0.4) of GRB afterglows. We compare the Gamma-ray 
burst Optical and Near-infrared Detector (GROND) GRB afterglow sample - one of the largest and most homogeneous samples of 
GRB fields - with galaxies and AGN found in the COSMOS-30 photometric catalog. We find no significant signal of anomalous 
clustering of galaxies at an estimated median redshift of z ~ 0.3 around GRB sight-lines, down to K AB < 19.3. This result is contrary 
to the expectations from the Mgll excess derived from GRB afterglow spectroscopy, although many confirmed galaxy counterparts 
to Mgll absorbers may be too faint to detect in our sample - especially those at z > 1 . We note that the addition of higher sensitivity 
Spitzer IRAC or HST WFC3 data for even a subset of our sample would increase this survey's depth by several orders of magnitude, 
simultaneously increasing statistics and enabling the investigation of a much larger redshift space. 

Key words. Gamma-ray burst: general - Galaxies: statistics 
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1. Introduction 

Gamma-ray bursts (GRBs) have proven to be powerful tools for 
studying the high redshift universe. Their afterglows pinpoint 
star forming galaxies that would otherwise be exceedingly diffi- 
cult to discover. Using bright transients such as GRBs to study 
the Universe introduces completely different selection criteria 
from standard surveying techniques. Long duration GRB (in 
this work, the term "GRBs" implies long-duration GRBs, un- 
less explicitly stated otherwise) afterglows have been used to 
study the star formation and metallicity evolution of the Universe 
(see e.g. Savaglio et al. 2009; Jakobsson et al. 2005; Prochaska 
et al. 2006; Salvaterra et al. 2012). The host galaxies of z > 5 
GRB afterglows have been studied by Tanvir et al. (2012); Basa 
et al. (2012). The upper limits on the non-detections that the au- 
thors derive suggest that the galaxy luminosity function evolves 
rapidly at these higher redshifts. In general, GRBs seem to be 
clear tracers of star formation (see e.g. Le Floc'h et al. 2006; 
Butler et al. 2010), however it is unclear the extent to which 
selection effects against more dusty and metal-rich galaxies may 
affect surveys(Kriihler et al. 201 1 ; Elliott et al. 2012; Perley et al. 
2013). 

In addition to studies of the GRB host galaxies, GRB after- 
glows themselves offer a brief and bright glimpse into the high 
redshift universe. Intervening absorption line systems can be 
detected in the same manner employed during the past decades 
with quasar spectroscopy. One such absorption line system 



is Mgll, which is easily detected in moderate S/N spectra at 
Arest ~ 2800 A. Because of its strong absorption and ease of 
identification owing to the fact that it is an absorption doublet, 
Mgll has been used extensively as a tracer of galaxies, galac- 
tic outflows, and chemical evolution. Mgll itself is coincident 
with a wide range of neutral hydrogen column densities from 
N m ~ 10 16 - 10 22 cnr 2 (Churchill et al. 2000), though there is 
strong evidence that Mgll equivalent width (EW) is directly cor- 
related with higher hydrogen columns and thus smaller impact 
parameters to the bulk star forming region of galaxies (Steidel 
et al. 1995; Bouche et al. 2006; Kacprzak & Churchill 2011; 
Bordoloi et al. 201 1). Indeed, Bordoloi et al. (2012) showed that 
Mgll absorbers follow a bi-modal spatial distribution, wherein 
at impact parameters smaller than d ~ 40 Kpc, Mgll is associ- 
ated with cool star forming regions, and at d > 40 Kpc Mgll 
tends to be uniformly distributed around galaxies, perhaps sug- 
gesting large scale outflows. There is no evidence to suggest 
that this behavior evolves with redshift, at least in the interval 
0.4 < z < 2. It is, however worth noting that Matejek et al. 
(2013) have studied the 2 < z < 6 regime using infrared spectra 
of quasars. The authors find that the association of Mgll with 
damped Lyman-a systems (DLAs) strengthens towards higher 
redshifts, even though the overall taxonomy of these absorbers 
as defined in Churchill et al. (2000) does not evolve. There is ev- 
idence that the redshift evolution in the number density of Mgll 
absorbers follows the cosmic star formation rate (Zhu & Me- 
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nard 2012; Matejek & Simcoe 2012), though the effects of ob- 
servational bias may still play an important role (Lopez & Chen 
2012). 

Prochter et al. (2006) compared the incidence of strong Mgll 
absorbers in quasar sight-lines with those found in GRB after- 
glows. In this context, "strong" is defined as absorption sys- 
tems whose rest frame equivalent width of the /12796 feature is 
Wo > 1.0 A. The authors found a factor of ~ 4 excess in the 
number density of these strong intervening 0.4 < z < 2.0 Mgll 
systems in GRB sight-lines relative to those found in quasars 
sight-lines. More recently, Vergani et al. (2009) have confirmed 
a factor of ~ 2 excess with a much larger sample. Besides in- 
creasing the statistical significance of the Mgll discrepancy, they 
found that the properties of weak 0.3 < Wo < 1 .0 A intervening 
Mgll systems are statistically identical between the two types of 
sight-lines. Interestingly, the abundance of CIV - a higher ion- 
ization absorption line system - also does not show an overden- 
sity (Sudilovsky et al. 2007; Tejos et al. 2007). Many authors 
have examined possible solutions to the observed discrepancy, 
and agree that a satisfactory explanation does not yet exist (Frank 
et al. 2007; Sudilovsky et al. 2009; Cucchiara et al. 2009; Ver- 
gani et al. 2009; Kann et al. 2010; Wyithe et al. 201 1; Rapoport 
et al. 2012). A comprehensive overview of these possible solu- 
tions is given in Porciani et al. (2007). It is especially intriguing 
that the amplitude of the discrepancy seemingly depends on the 
resolution of the afterglow spectra suggesting either that rapid 
response high resolution spectroscopy may introduce a not yet 
understood observational bias in any sample study of GRB af- 
terglows, or that the Mgll discrepancy is a statistical fluke (Cuc- 
chiara et al. 2012). 

Objects near the line of sight to a GRB afterglow are un- 
affected by light from the event after the afterglow has faded, 
allowing detailed follow-up of the field to extremely high an- 
gular resolutions. The study of sources near quasars are re- 
liant on either extremely high resolution imaging on intrinsically 
dim quasars, Lyman-a imaging, or very low redshift quasars. 
Galaxies that give rise to absorption line systems in GRB after- 
glow spectra have been directly imaged and studied, although the 
number of robust associations is still extremely small. Chen et al. 
(2009) find that additional galaxies are found at very close angu- 
lar distances to GRB host galaxies whose afterglows exhibited 
strong Mgll absorption, though it is unclear what fraction were 
associated with the absorption line systems. Schulze et al. (2012) 
proposed galaxy counterpart candidates to absorbers, and addi- 
tionally found field galaxies with the same redshift as absorption 
line systems at distances of 130-161 kpc away from the sight- 
line. The frequency of field galaxies that are associated with 
Mgll absorption line systems is still unclear. However, Lopez 
et al. (2008) have shown that strong Mgll absorbers are strongly 
associated with galaxy clusters, and that the number density of 
these absorbers is much higher relative to field galaxies. In this 
work, we test for any anomalous signature of clustering or over- 
abundance of field galaxies in GRB and quasar lines of sight. 

If there is indeed a higher probability of detecting strong 
Mgll absorption line systems in GRB afterglow than in quasar 
sight-lines, one should also expect to detect more galaxies at 
close angular separations to GRB sight-lines. The two point cor- 
relation function is a powerful tool to determine if objects are 
clustered, and, if so, what their correlation lengths are. The an- 
gular two point correlation function is formally defined as 



dP = n[l + w(0)]d£l, (1) 
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where dP is the probability of finding an object within a solid 
angle c/Q at an angular distance 8 (Peebles 1980). In general, a 
positive w implies some enhancement of object-object grouping 
above a uniform random distribution, while a negative w implies 
some avoidance. The two point correlation function has been 
used to estimate the clustering properties of galaxies, quasars, 
and GRBs (see e.g. Groth & Peebles 1977; Ross et al. 2009; 
Brainerd et al. 1995). 

In this work, we measure the angular two-point correlation 
function between GRB afterglows and field galaxies, and com- 
pare this quantity with galaxy-galaxy and AGN-galaxy correla- 
tions. In §2, we present the sample that we analyze in §3. We 
present the results of the analysis in §4, and discuss these results 
in §5. 

2. The sample 

The Gamma-ray burst Optical and Near-infrared Detector 
(GROND) is a simultaneous 7-channel imager mounted on the 
MPG/ESO 2.2m telescope at La Silla, Chile (Greiner et al. 
2008). The four optical channels are nearly identical to the SDSS 
g'r'i'z' filters, while the infrared channels are effectively equiv- 
alent to the 2MASS JHK filters. While the detector system was 
built to quickly measure photometric redshifts of GRBs via the 
Lyman-a dropout technique, the multi-wavelength photometry 
provided by GROND can potentially reveal a wealth of physical 
characteristics of any source (Greiner et al. 2008). 

Our sample of GRB fields is constructed exclusively of 
GROND observations taken between 2007-2012. This fact en- 
sures a high degree of homogeneity, as no correction for cross- 
instrumental calibrations are required. Furthermore, the simulta- 
neous nature of GROND observations ensures that data on a per- 
field basis have not varied due to intrinsic variability or weather 
conditions. Since only ephemeris, hardware failure and weather 
losses prevent GROND follow-up observations, the GROND 
sample of GRB afterglows has an exceptional success rate of 
~ 90% detection of long-duration GRBs when follow-up is pos- 
sible within the first few hours after the trigger (Greiner et al. 
2011; Kriihler et al. 2011). We include in our sample all fields 
which 1) have a XRT localization of the burst, 2) Galactic lat- 
itude \b\ > 10°, 3) at least 15° from the Galactic center, and 
4) have been observed to a 3 cr background limiting magnitude 
of Kab ~ 20.4. This limiting magnitude corresponds roughly 
to an integration time of 1500 and 1200 seconds in optical and 
NIR, respectively, among four telescope dither positions. For 
reference, typical 3<x AB limiting magnitudes in the optical are 
at least 3 mag deeper than K, while J and H are 1.0 and 0.5 
mag deeper, respectively. Our final sample consists of 73 GRB 
fields (see Table 1 and Fig. 1). Since most of the bursts have no 
measured redshift, we assume that our bursts follows the same 
redshift distribution as in the TOUGH sample (Jakobsson et al. 
2012). 

3. Procedure 

3.1. Creating the source catalog 

To accurately measure the optimal extraction region and measure 
the apparent extent of sources from GROND multi-band images, 
we co-add all 7 images to create a single detection image for ev- 
ery field. Since the exposure time differs between the NIR and 
optical channels, we scale and weigh the images to a common 
effective exposure. SExtractor (Bertin & Arnouts 1996) is run 
on the detection image using relatively low source extraction 
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Fig. 1. Positions in equatorial coordinates of the 73 GROND GRB 
afterglows in the sample. Bursts within 10° of the Galactic plane are 
excluded. 

thresholds such that 5-10% of the sources are likely spurious 
detections, but still no bright sources are incorrectly de-blended 
or otherwise split into multiple sources. The detection image 
source catalog is correlated with each individual band's source 
catalog, and only those sources whose 1) positions in each cata- 
log are consistent to within 0'.'5 of each other, and 2) are detected 
in at least 6 out of the 7 bands (one of which is required to be a 
K-band detection), and 3) have an error <tk < 0.10 mag are used 
in the final source catalog for that field. The latter criteria limits 
the overall sample to roughly cr K < 19.3, with ~ 10% of sources 
dimmer than this. The incompleteness down to K ~ 20.0 does 
not introduce any significant bias in the final results. This cross- 
correlation ensures that spurious detections are removed, and 
that each source has multi-wavelength detections. The require- 
ment that objects be strongly detected in K ensures that they can 
be reliably seperated between galaxies and stars. 

Creating a source catalog based on a detection image in the 
aforementioned manner minimizes the probability that flux from 
extended objects is missed due to varying spectral properties as a 
function of position. That is, by determining a suitable aperture 
from the co-added image, one is guaranteed that the entirety of 
flux in an extended source is included therein if that same aper- 
ture is then used in each individual image, as long as that source 
is strongly detected in each band. Additionally, if the source is 
reasonably isolated in the image, this method avoids the com- 
plications introduced by PSF matching and galaxy fitting while 
still providing adequate photometry. We therefore perform aper- 
ture photometry with SExtractor in dual-mode using the detec- 
tion image as the template for each band. 

The SExtractor magnitudes are calibrated based on time- 
tabulated instrumental zeropoints in the case of the optical chan- 
nels and 2MASS field stars in the case of the NIR. GROND ze- 
ropoints are calibrated on average once every three months based 
on SDSS standard fields, and immediately after technical work 
is performed on the optical system. The RMS of the optical ze- 
ropoints as a function of time is on the order of 0.05 mag. We 
expect a similar spread due to varying atmospheric conditions, 
as the range of seeing in our sample is between 1.0 and 1.6 in r'. 
In 18 fields (24%) co-incident with the SDSS, magnitudes were 
calibrated directly against SDSS field stars. 

Next, we fit each object's spectral energy distribution (SED) 
with LePHARE (Arnouts et al. 1999; Ilbert et al. 2006). LeP- 
HARE is a spectral template fitting tool based on x 2 minimiza- 
tion. We determine the best fit galaxy and stellar templates for 
each source, using the COSMOS galaxy templates, with emis- 
sion lines and prescribed reddening and parameters therein (Il- 
bert et al. 2009). More specifically, we adopted the Prevot (Pre- 
vot et al. 1984) extinction law for late type templates and the 
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Fig. 2. A" A b magnitude distribution of the 3068 field galaxies identified 
in the sample. The steep drop-off after K ~ 19.3 is due to the criteria 
that o~k < 0. 10 for each galaxy, and that is it detected in five other 
bands. This represents a decline in completeness of the sample. The 
typical 3<x limiting magnitude for an individual image as determined 
from the sky background is K A % ~ 20.4. 



Calzetti law (Calzetti et al. 1994) and two modifications thereof 
(Ilbert et al. 2009) for the SB templates. The values of extinc- 
tion range from 0.0 < E(B - V) < 0.5 in steps of 0.1. Early type 
templates are not corrected for extinction as there are not well 
tested empirical models to describe the dust distribution for these 
galaxies. In the case of galaxy templates, redshift is constrained 
to z < 2 with a step size Az = 0.01. The former constraint 
removes some degeneracy between high and low redshift tem- 
plate fits, and is justified since the population of galaxies with 
an apparent Kab < 20.0 is negligibly small at z > 1.5. We as- 
sume that the contamination from AGN-dominated sources in 
our galaxy catalog is negligible (< 5%, as discussed later in the 
paper). We separate galaxies and stars based on their best fit tem- 
plate and shape parameters. Any sources with a full-width half 
maximum (FWHM) or ellipticity as measured with SExtractor 
in 5<x excess of the stellar values for that field are automatically 
classified as galaxies. These stellar values are determined by 
computing the average and standard deviation of the lowest 20% 
FWHM and ellipticity distributions, under the assumption that 
those "bottom 20%" of sources are stars. In practice, this usually 
means the FWHM and ellipticity of 5-10 stars are averaged for 
a given field. Besides this shape criterion, we categorize sources 
with 1 .5 x xlm > Afgaia X y as galaxies, where x 2 corresponds to 
the best fit template for that class of source. This latter crite- 
rion has been applied to the COSMOS field with great success 
(Salvato et al. 2009). The extinction law and extinction values 
are to some extent free parameters in the fitting. Therefore, a 
certain amount of degeneracy in redshift-color-extinction space 
is expected. However, we do not expect this degeneracy to af- 
fect star/galaxy classification. This is due to the consideration of 
morphological information (i.e. extended vs point-like) and the 
implied high redshift limit imposed by the K band magnitude 
limit. 

We detect 3068 galaxies and 1368 stars that match our crite- 
ria. The K-band magnitude distribution of galaxies is presented 
in Fig. 2. 
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Fig. 3. Color-color diagrams for stellar templates (open circles) and sources identified as stars (black triangles and blue squares) in the GROND 
sample. Offsets of +0.2, +0.08, +0.08, and +0.1 mag to the q' , t ', i', and z' zeropoints are applied, except in the case of direct SDSS calibration. 
The sources represented by black triangles are SDSS calibrated, while the blue squares represent a zeropoint calibration. 



3.2. Verifying the source catalog 



To determine the accuracy of our galaxy-star separation, we 
compare a catalog of sources derived from nine GROND sub- 
fields of the COSMOS field with the COSMOS 30-band photo- 
metric catalog (COSMOS-30; Ilbert et al. 2009). COSMOS-30 
has much larger spectral coverage and higher sensitivity than the 
GROND observations, making it an ideal standard with which 
to compare. We employ COSMOS-30 as both a measure of 
the quality of galaxy identification from GROND and to calcu- 
late the AGN-galaxy two point correlation function. From the 
nine adjacent sub-fields observed by GROND (center pointing 
at R.A.=150.0862 deg, Dec.=2.3745 deg), we identify 317 total 
candidate galaxies. Nine of these candidates are in a masked area 
in COSMOS-30, and a further nine are categorized as stars in 
COSMOS-30. This corresponds to 2.9% of galaxies that are mis- 
classified using our method, assuming COSMOS-30 has 100% 
accuracy in classification. Besides the masked objects, no candi- 
dates were detected in GROND images without a corresponding 
source in COSMOS-30. The lack of spurious detections is likely 
a result from our stringent detection criteria discussed in §3.1, 
namely that an object must have independent detections in at 
least six filters. 



We evaluate our photometric calibration by comparing the 
colors of sources that we classify as stars to the stellar templates 
of Pickles (1998) and Bohlin et al. (1995). In Fig. 3, we present 
color-color diagrams of our sources compared with those of stel- 
lar templates. On average, offsets of +0.2, +0.08, +0.08, and 
+0.1 mag to the g', f , i', and z' zeropoints, respectively, are re- 
quired to match the colors of the stellar templates. In the case 
of direct SDSS calibration, no offset is needed nor applied. Af- 
ter applying these offsets, the colors of stars in our catalog are 
on average not offset from those of stellar templates, implying a 
reliable photometric calibration. The somewhat higher scatter of 
our objects relative to the templates is due to a combination of 
uncertain dust extinction towards foreground stars, binary sys- 
tems, and photometric uncertainty. 

3.3. Measuring the two-point correlation 

After the source catalog is categorized, calculating a two-point 
GRB-galaxy correlation function is possible. Since our sample 
consists of isolated fields with no overlap, we must calculate dis- 
tance pair distributions on a field by field basis, later combining 
them into a global correlation function. This technique limits 
the angular sizes we can examine, and furthermore introduces 
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systematic errors when the size scales approaches the size of in- 
dividual images. However, since we are specifically examining 
clustering at small scales, GROND's 5.4' x 5.4' images 1 have 
sufficient spatial coverage. We limit our our analysis to angular 
distances < 120". This angular distance corresponds to ~ 600 
kpc at z — 0.4, assuming a ACDM concordance cosmology with 
h = 0.7 and Q m = 0.27 (Wright 2006). 

Many of the afterglows in our sample have no measured red- 
shift. We assume that our sample of GRB afterglows has an un- 
derlying redshift distribution that is consistent with the literature 
(Jakobsson et al. 2012). More specifically, the GRBs should be 
background to the galaxies in the field. This latter assumption is 
likely justified, since the mean redshift of GRBs is z ~ 2 and the 
mean redshift of galaxies with K < 19.3 is z ~ 0.4, as noted in 
§2. 

To measure the two-point correlation, we start by calculating 
the angular distance between the GRB optical afterglow posi- 
tions (or the XRT afterglow position circle, in the case that no 
optical afterglow was detected) and each galaxy. This provides a 
distribution of arclengths hod for a given field. Next, we assign 
a new position to each galaxy within that field based on values 
picked from a uniform random distribution. Any galaxies that 
are within two stellar FWHM of a star or fall into an object mask 
are re-assigned a position until they do not. This ensures that 
the observability of the isotropically distributed control sample 
shares the same observational properties as the actual galaxies. 
We calculate the distances between each GRB-random galaxy 
n DR and the random sample's mutual separation n RR . This pro- 
cess is repeated and averaged over a total of 10 3 times for each 
field. We then calculate the angular correlation function w and 
its variance cr using the Landy & Szalay (Landy & Szalay 1993) 
estimator 



DD - 2DR + RR 

W = RR ' (2) 

oi= 1 -^, (3) 

n DD 

where DD, DR, and RR are the normalized frequency dis- 
tributions of tiDD, nDR, and n RR , respectively. We perform 
the same procedure to measure w ga i axv - ga i axy , w AC N-gaiaxy, and 
^random-random- Though it is trivial to show analytically that w 
is zero in the case of isotropically distributed data, we explicitly 
perform the calculation using three randomly chosen coordinates 
per field as a verification of our method. 

We compute the correlation functions on both the GROND 
and COSMOS-30 sample. To approximate the same selec- 
tion criteria and completeness of our GROND sample, we re- 
move from the COSMOS-30 catalog any sources dimmer than 
Kab = 19.3 as well as those sources whose photometry may be 
unreliable due to saturation or being within an object mask (No 
flags in any optical filter). This leaves 4481 galaxies and 2273 
stars. The photometric redshift distribution of these COSMOS- 
30 galaxies is presented in Fig. 4. The mean and median redshift 
of this sub-sample is 0.387 and 0.361, respectively. For comput- 
ing wxgn -galaxy, we ensure that only background X-ray selected 
sources are included in the calculation by imposing a z > 1.0 
redshift criterion (Salvato et al. 2011). This criterion limits the 
catalog to 980 sources. 



1 5.4' x 5.4' in g'r'i'z', 10.3' x 10.3' in JHK 
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Fig. 4. Photometric redshift distribution of the 448 1 galaxies in 
COSMOS-30 that have K AB < 19.3. The mean and median redshifts 
of the distribution are 0.387 and 0361, respectively. 

4. Results 

In Fig. 5, we present the two point GRB-galaxy, galaxy-galaxy, 
and random-random angular correlation functions computed 
from the GRB sample. The two point correlation functions are 
self-consistent for all angular distance scales, with the exception 
of the (as expected) null random-random function. A possible 
contribution to the first bin of the GRB-galaxy wgrb-gal is likely 
due to the detection of the host galaxy of the GRB, although the 
vast majority of known GRB hosts are fainter than our magni- 
tude limit. This corresponds to 3" in our sample, and it is en- 
tirely possibly to find GRB host galaxies at up to these distances 
from the afterglow position and larger. Furthermore, these are 
also the angular distance scales within which Chen et al. (2009) 
find sources nearby sight-lines to GRBs with spectroscopically 
identified intervening Mgll systems, albeit at very faint magni- 
tudes. 

We also find clustering of galaxies within small angular dis- 
tances of each other as seen in the galaxy-galaxy correlation 
wgal-gal- This is consistent with expectations from multiple 
galaxies occupying the same massive halo (Berlind & Weinberg 
2002; Lee et al. 2006) and naturally from galaxy clusters. This 
signal is seen more clearly in the correlation function as mea- 
sured from the COSMOS-30 sample (Fig. 6), where we include 
only sources with K AB < 19.3 to approximate our sample's com- 
pleteness. For the COSMOS-30 sample, the best fit power-law to 
the first 14 bins (0 < 6 < 42") of the form (6/6qT s yields a cor- 
relation length # = 2.6 + 0.4 arcseconds and slope 6 = 0.8 ± 0.2 
(statistical error only). This same power-law is presented in all 
two-point correlation functions for comparison. 

The AGN-galaxy two-point correlation function w ACN _ CAL is 
also presented in Fig. 5. 

5. Discussion 

A wcRB-gai that is consistent with w C al-gal and w AG n-cal im- 
plies that there is no excess or unusual clustering of galaxies 
around GRB events, contrary to the expectation from observa- 
tions of Mgll absorption line systems. Though our data support 
this null hypothesis, one must consider that the limiting mag- 
nitude of our sample is relatively shallow in comparison to the 
expected magnitudes of galaxies hosting Mgll absorbers at the 
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Fig. 5. (Top to bottom) GRB-galaxy, galaxy-galaxy, and random- 
random two point correlation functions. The solid black line represents 
the best fit power-law to the COSMOS-30 galaxy-galaxy two-point cor- 
relation function, as described in Fig. 6 and in §4. The random-random 
two-point correlation function is, as expected, consistent with zero. 



redshifts found in GRB afterglow spectra. Indeed, the typical 
magnitudes of candidate absorption line counterparts detected 
by Schulze et al. (2012) have K AB > 22, and Chen et al. (2009) 
find that the objects at angular distances of sources 1" - 3" to 
GRB hosts with spectroscopically detected Mgll typically have 
H AB >26. 



Although we expect that an overdensity of Mgll absorbers 
corresponds to an overdensity of galaxies around the line of 
sight, we are unfortunately limited to the very brightest and clos- 
est objects in our survey. A qualitative measure of how much 
low number statistics affect w is demonstrated by comparing our 
galaxy-galaxy correlation with that calculated from COSMOS- 
30 (Figs. 5-6). The correlation as measured by 3068 galaxies in 
GROND fields only shows a hint of clustering due to halo occu- 
pation of multiple galaxies and galaxy clusters, as discussed in 
§ 4. The signal is much more clear from the 448 1 galaxies with 
K AB < 19.3 that we extract from COSMOS-30. 

Clustering analysis of GRB afterglows and/or hosts have 
been performed by Wang & Wei (2010); Bornancini et al. 
(2004). Wang & Wei (2010) find no evidence of clustering be- 
tween SWIFT XRT afterglows and ROSAT selected galaxy clus- 
ters. The purpose of our survey is not to examine large scale 
anisotropics, but rather the small scale clustering implied by the 
Mgll excess. Bornancini et al. (2004) examined whether GRB 
host galaxies tend to reside in high density environments by 
calculating the GRB-galaxy two point correlation for 19 GRB 
hosts. They concluded that GRB host galaxies likely do not oc- 
cur in over-dense areas, although local cosmic variance could 
still significantly affect these results due to the small sample size. 
Furthermore, choosing hosts over afterglows introduces an ad- 
ditional bias against intrinsically dim or reddened hosts. Con- 
structing a sample based on afterglow positions as we have done 
in this work removes this bias and increases the sample size sig- 
nificantly. It is furthermore justified, since there is no a priori 
reason to link foreground absorption in the afterglow with the 
detection or non-detection of a host galaxy. 

The largest limitation to our survey stems from the fact 
that there is significant degeneracy between galactic/galactic and 
galactic/stellar templates in GROND's optical regime for ob- 
jects at z < 2. For this reason, including NIR data is cru- 
cial to distinguish point-like galaxies from stars. The sensitiv- 
ity of GROND's K channel is typically three magnitudes shal- 
lower than its optical channels. Thus, we are forced to exclude 
a large number of sources from our survey simply due to insuf- 
ficient sensitivity. Including Spitzer/IRAC or HST/WFC3 data 
for a significant fraction of fields observed by GROND would 
remove the largest constraint to our survey by providing much 
deeper sensitivity limits and/or wavelength coverage, simulta- 
neously increasing statistics and enabling the investigation of a 
much larger redshift space. Such surveys are currently under- 
way and are scheduled to be completed in the next years (Levan 
2009b,a; Perley et al. 2012). Combining these high quality NIR 
data with GROND data would also provide for accurate pho- 
tometric redshifts, enabling a measure of the two point spatial 
correlation function. 
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6 (arcsec) 6 (arcsec) 



Fig. 6. (Left to right) AGN -galaxy and galaxy-galaxy two-point correlation functions computed from the COSMOS photometric catalog. To 
approximate the GROND sample selection and completeness, only sources with K AB < 19.3 are included. The best fit power-law of the first 14 
bins (0 < 6 < 42") of the galaxy-galaxy two-point correlation function is represented by the solid black line. The same fit is used in both figures 
for comparison. 
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Table 1. Identification and positions of the GRBs in the final sample. 



GRB 


R.A. J2000 


Dec. J2000 


071112C 


02: 


:36 


:50 


.93 


+28:22:16.68 


080212 


15 


:24 


:35 


.42 


-22:44:29.70 


080330 


11 


:17 


:04 


.51 


+30:37:23.48 


080408 


07: 


:38 


:39 


.59 


+33:18:14.90 


0804 13B 


21: 


44 


:34 


.67 


-19:58:52.40 


080413 


19: 


:09 


:11 


.74 


-27:40:40.30 


0805 14B 


21: 


:31 


:22 


.69 


+00:42:28.90 


080520 


18: 


40 


46 


.30 


-54:59:31.00 


080523 


01 


:23 


:11 


.70 


-64:01:51.50 


080605 


17: 


:28 


:30 


.05 


+04:00:55.97 


080710 


00: 


:33 


:05 


.67 


+ 19:30:05.69 


080916 


22: 


:25: 


:06 


.20 


-57:01:22.90 


080928 


06 


:20 


:16 


.82 


-55:11:59.30 


081008 


18: 


:39: 


49 


.89 


-57:25:52.80 


081109 


22 


:03 


:09 


.86 


-54:42:41.04 


081118 


05: 


:30 


:22 


.18 


-43:18:05.30 


081121 


05: 


:57 


:06 


.08 


-60:36:09.94 


081228 


02: 


:37 


:50 


.89 


+30:51:09.10 


081221 


01 


:03 


:10 


.19 


-24:32:52.20 


081230 


02 


:29 


:19 


.53 


-25:08:51.72 


090102 


08: 


:32 


:38 


.10 


+33:11:45.30 


090123 


00: 


:27 


:08 


.74 


-23:30:04.00 


090205 


14: 


43 


:38 


.65 


-27:51:10.70 


090323 


12: 


42 


:50 


.29 


+ 17:03:11.60 


09040 IB 


06 


:20 


:21 


.10 


-08:58:19.35 


090423 


09 


:55 


:33 


.29 


+ 18:08:58.00 


090424 


12: 


:38 


:05 


.09 


+ 16:50:14.75 


090426 


12: 


:36 


:18 


.04 


+32:59:09.24 


090509 


16 


:05 


:39 


.01 


-28:23:59.64 


090516 


09: 


:13 


:02 


.59 


-11:51:14.90 


090519 


09 


:29 


:07 


.00 


+00:10:49.10 


090530 


11: 


:57 


40 


.51 


+26:35:38.40 


090812 


23: 


:32 


48 


.54 


-10:36:17.60 


090814 


15: 


:58 


:26 


.35 


+25:37:52.42 


090902B 


17: 


:39 


45 


.41 


+27:19:27.10 


090926B 


03: 


:05 


:13 


.94 


-39:00:22.20 


091018 


02: 


:08 


44 


.61 


-57:32:53.70 


091029 


04: 


:00 


42 


.60 


-55:57:19.84 


091109 


20: 


:37 


:01 


.80 


-44:09:29.60 


100219A 


10: 


:16 


48 


.50 


-12:34:00.50 


100414A 


12: 


48 


:26 


.93 


+08:41:34.40 


100518A 


20: 


:19 


:09 


.33 


-24:33:16.56 


100621A 


21 


:01 


:13 


.11 


-51:06:22.46 


100902A 


03: 


:14: 


:30 


.96 


+30:58:45.23 


101023A 


21: 


:11 


:51 


.23 


-65:23:15.61 


101219B 


00: 


48 


:55 


.34 


-34:33:59.26 


110128A 


12: 


:55 


:35 


.10 


+28:03:54.10 


1 10206A 


06 


:09 


:20 


.04 


-58:48:24.91 


110312A 


10 


:29 


:55 


.47 


-05:15:45.20 


1 10407A 


12 


:24 


:07 


.49 


+ 15:42:42.16 


110709B 


10 


:58 


:37 


.11 


-23:27:16.76 


110721A 


22: 


:14 


:38 


.19 


-38:35:35.70 


110818A 


21: 


:09 


:20 


.89 


-63:58:51.80 


111008A 


04: 


:01 


48 


.22 


-32:42:34.09 


111107A 


08: 


:37 


:54 


.65 


-66:31:12.40 


111129A 


20: 


:29 


44 


.14 


-52:42:46.48 


11 1209 A 


00 


:57 


:22 


.70 


-46:48:05.00 


111211A 


10: 


:12 


:21 


.70 


+ 11:12:30.00 


111212A 


20: 


:41 


43 


.52 


-68:36:45.00 


111228A 


10 


:00 


:16 


.01 


+ 18:17:51.80 


11 1229 A 


05 


:05 


:08 


.84 


-84:42:38.70 
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1201 19A 


08:00 


:06.94 


-09:04:53.83 


12021 1A 


05:51: 


:00.89 


-24:46:30.79 


120302A 


08:09: 


:35.54 


+29:37:41.05 


12031 1A 


18:12: 


:22.16 


+ 14:17:46.30 


120320A 


14:10 


:04.30 


+08:41:47.26 


120404A 


15:40. 


:02.29 


+ 12:53:06.29 


120422A 


09:07 


:38.42 


+ 14:01:07.36 


121024A 


04:41 


:53.30 


-12:17:26.48 


121027 A 


04:14 


:23.45 


-58:49:47.17 


121217A 


10:14 


:50.51 


-62:21:03.28 


121229 A 


12:40: 


:24.29 


-50:35:39.48 
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